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Abstract 

We examine the stability of the "coast" motion of fish, that is to say, the motion of a neutraUy buoyant 
fish at constant speed in a straight hne. The forces and moments acting on the fish body are thus perfectly 
balanced. The fish motion is said to be unstable if a perturbation in the conditions surrounding the fish 
results in forces and moments that tend to increase the perturbation and it is stable if these emerging 
forces tend to reduce the perturbation and return the fish to its original state. Stability may be achieved 
actively or passively. Active stabilization requires neurological control that activates musculo-skeletal 
components to compensate for the external perturbations acting against stability. Passive stabilization 
on the other hand requires no energy input by the fish and is dependent upon the fish morphology, i.e. 
geometry and elastic properties. In this paper, we use a deformable body consisting of an articulated 
body equipped with torsional springs at its hinge joints and submerged in an unbounded perfect fluid 
as a simple model to study passive stability as a function of the body geometry and spring stiffness. 
We show that for given body dimensions, the spring elasticity, when properly chosen, leads to passive 
stabilization of the (otherwise unstable) coast motion. 



1 Introduction 

We analyze the passive stability of a class of elastic bodies moving in a perfect fluid. Our primary motivation 
is to examine the role of body elasticity as a mechanism for passively stabilizing the coast motion of fish. 
Fish seem to alternate between actively swimming ( "burst" ) and passively responding ( "coast" ) to the 
surrounding fiuid - a cycle known as the "burst and coast" cycle, see [26]. It is of course difficult to gauge 
whether the coast motion is entirely passive, that is, whether or not the fish employs active control for 
stabilization purposes; see, for example, '27]. Yet, experimental evidence seems to suggest the existence 
of some mechanisms for passive stabilization in fish motion. For example, [3] reported upstream motion 
of anesthetized fish (undeniably passive) in vortical flows and mentioned that the elastic properties of the 
fish body are essential to achieve this motion. Indeed, the motion of elongated rigid bodies in unsteady 
flows is known to be passively unstable, see, for example, |151 IHl [20]. Body elasticity may not be the only 
mechanism responsible for passively stabilizing such motion. Other physical mechanisms such as the effects 
of the fish wake or the boundary layer effects may also be at play. In this paper, we focus solely on the role 
of body elasticity in the passive stability of the coast motion. This is of course a mathematical idealization 
while, in reality, the effects of body elasticity, viscosity and vortex shedding are all coupled and operate 
simultaneously. This mathematical idealization is aimed to untangle the effects of body elasticity on the 
stability of motion and is viewed as a first step in constructing models with higher degrees of complexity 
and fidelity to the physical system. 

We model the fish as an articulated body formed of {n + 1) rigid links, say, identical ellipses of major- 
and minor-axes a and 6, respectively. The links are connected via n hinge joints with torsional springs at 
the joints (with stiffness fc^,? = l,...,n) to emulate the elasticity of the fish body. The articulated body 
is submerged in an unbounded volume of incompressible, inviscid fluid at rest at infinity. The dynamics 
of this body-fluid system admits a family of relative equilibrium solutions where the body, in its elongated 
configuration, is steadily moving along its major axis of symmetry, see figure[T] We refer to these translational 
motions as the coast motion. We examine the effects of the springs' elasticity on the stability of the coast 
motion. We focus on passive stability, that is to say, the stability of motion when the deformable body is 
subject to initial conditions only and is passively responding to the surrounding fluid. This is in contrast to 
a deformable body that is actively controlling its shape to swim (see ^) or to stabilize its motion (see [27j). 
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(a) unstable for all b/a < I (b) passively stable for a range of fe/a < 1 and k values 




(c) passively stable for a range of 6/a < 1 and fci, k2 values 



Fi gure ll Family of relative equilibria for articulated body moving in an inviscid, incompressible fluid: (a) rigid ellipse is 
unstable for b < a; (b) and (c) two-link and three-link bodies may be stable for proper choice of spring stiffness versus body 
dimensions. 

While the motion of a single rigid body, see figure [ija) is known to be passively unstable, we find that, for 
the elastic body, there exists a range of parameter values - spring stiffness versus body geometry - for which 
the coast motion is passively stable. We examine the regions of stability for a two-link and a three-link body. 
We find that the stable region of the two-link body is characterized by a bending mode of deformations. For 
the three-link body, there exists two stable regions, one is characterized by a bending mode of deformation 
and the second by a traveling-wave mode of deformation. 

This problem is somewhat reminiscent to that of the fluttering flag instability. The latter has been the 
focus of a large number of experimental ([5T]), computational ([T1[TH]), and theoretical investigations. 
Of course, the flag is held stationary at one end in an incoming uniform flow fleld whereas the fish is freely 
moving. This constraint does alter the forces applied on the flag and, in turn, may significantly influence the 
stability results. For example, [22] and [12] showed that the stability of the freely moving Foppl equilibria is 
distinct from the stability of the stationary equilibria in an incoming uniform flow fleld. Further, the flag is 
assumed to be heavy whereas the fish is neutrally-buoyant. The hydrodynamic forces acting on the flag are 
thus dominated by vortex shedding rather than the added mass effect. All is to say that, while the fish and 
flag are somewhat conjugates, the fluttering flag instability (typically analyzed as a function of the flag's 
rigidity versus weight) cannot be readily translated to the stability of the moving fish which we analyze in 
this work as a function of the body stiffness versus geometry. 

2 Problem Setting 

Consider an articulated body moving in an unbounded volume of incompressible, inviscid fluid at rest at 
infinity. The articulated body is composed of n + 1 identical ellipses (with semi-axes a and b) of uniform 
density p equal to that of the fluid. That is to say, the articulated body is neutrally-buoyant. The ellipses 
or links forming the body are connected via n massless and frictionless hinge joints placed at a distance I 
from the ellipses' centers along the major axes and equipped with torsional springs of constant stiffness ki, 

I = 1, . . . , 71. 

Let {ei 2,3} be an orthonormal inertial frame where {ei 2} spans the plane of motion. The configura- 
tion of the articulated body can then be described by 2{n + 1) position coordinates {xj,yj) of the ellipses' 
centers and {n + 1) orientation variables /3j (where j = 0, . . . ,n) defined as the angles between the ellipses' 
major axes of symmetry and the ei-axis, see figure [2j This set of 3(n -I- 1) coordinates is subject to 2n 
holonomic constraints due to the presence of the hinge joints connecting the ellipses. These constraints are 
implicitly accounted for when describing the body's configuration in terms of {xo,yo): Po: and n relative 
angles 6i = /3i+i — (3i { i = 1, . . . ,n) that denote the angles between every two consecutive links of the body. 
The coordinates {xo,yai M coincide with the inertial coordinates of the first ellipse and are called locomo- 
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tion coordinates while the relative angles 9i describe only the shape of the articulated body. For notational 
convenience, we introduce the (n + 3) -dimensional configuration vector q = {xo,yo, l3o,9i, . . . ,0„}^, where 
0"^ denotes the transpose and the n-dimensional shape vector = {6i, . . . ,0n}^- 

It is convenient for deriving the equations of motion in S|3]to write the velocity of the articulated body 
in terms of an orthonormal body-fixed frame {bi^2,3} (with = e^) whose origin is attached at (xq, yo) and 
such that Po is the angle between bi and ei. Let {uj, Vj) denote the components of the translational velocity 
vector Vj {j — 0, . . . ,n) of the center of each ellipse expressed in the body frame {bi 2,3}- In other words, let 
Vj = Ujhi + Vjh2- Also, let = fljh^ denote the angular velocity vector of each ellipse. One has Hq = $q 
and flj = $0 + ^21=1 f"^^' i = 1; ■ • ■ : Here, the dot notation () is used to denote derivative with respect to 
time t. We now construct a 3(n -I- l)-dimensional velocity vector ^ = {uq, vq, ^Iq, mi, wi, f2i, . . . , m„, ?;„, fln}^ . 
The velocity vector ^ can be related to the velocity vector q via a 3(n -I- 1) x (3 -I- n) transformation matrix 
M, function of /3q and only, namely, 

I = m- (1) 

For example, in the case of an articulated body made of two ellipses connected via one hinge joint, the 6x4 
transformation matrix M is 
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(2) 



This transformation matrix can be readily generalized to arbitrary but finite n. 



3 Equations of motion 

We derive the equations of motion for the fluid-mass system using Hamilton's least action principle for which 
the Lagrangian function is given by 

C = T -V (3) 

where T = Tb + Tp denotes the kinetic energy of the fluid-mass system, whereas V denotes the potential 
energy stored in the springs and can be written in matrix form as 

V=^&^K&, (4) 

where the stiffness matrix K is an n x n diagonal matrix with entries ki corresponding to the stiffness of the 
torsional springs at the hinge joints. 
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The kinetic energy Tg of the articulated body is given by 





Ts-^I^Mbodyl, (5) 

where Mbody isa3(n + l)x3(n + l) diagonal mass matrix with 3x3 block diagonal entries of the form 

\ 

m-ciiipso I , (6) 


where TOoiUpsc = pi^ab and JoUipsc — pT^ab{o? + b"^)/^ are the actual mass and moment of inertia of each 
ellipse. 

In potential flow, the kinetic energy of the fluid Tp can be written, using standard vector identities and 
techniques (see for example 111] and [19! ), ^ function of the configuration and velocity of the submerged 
body, 

Tf = i^^Maddcd^ . (7) 

The added mass matrix Madded is a 3(n + l) x 3(n+l) symmetric matrix that accounts for the presence of the 
fluid. Differently said, in potential flow, the hydrodynamic forces and moments acting on the surfaces of the 
accelerating bodies can be completely accounted for by 'adding mass' to the submerged bodies (see ^7 for a 
comprehensive review). The entries of the added mass matrix Madded depend on the shape of the articulated 
body which may be changing in time as dictated by the evolution of the shape variable 0. We compute the 
entries of Madded (©) numerically using a panel method, as discussed in |111 113j . It is worth noting that, 
although the added mass effect captures all the hydrodynamic forces in the potential flow environment, this 
is not the case in real viscous fluids where one needs to account for additional drag forces and forces induced 
by vortex shedding and vortex wakes as discussed in ^jsj 

Adding ^ and ([7| and substituting ([T]) into the resulting expression, the total kinetic energy T — 
Tg + Tp of the fluid-mass system can be written as 

T=iq^Mq, M -M^ (Mbody + Madded (0))M, (8) 

where M is an (n + 3) x {n + 3) mass matrix as a function of /3o and only since the transformation matrix 
M depends only on /Sq and 0. The Lagrangian function can then be expressed as 

^(q,q) = ^q^Mq-i0^K0. (9) 

In the absence of external forces and moments acting on the fluid-body system, the equations of motion are 
given by the (n + 3) Euler-Lagrange equations. 

This class of problems is energy-preserving, hence the total energy E = -q-^Mq+ -0-^K0 is conserved. 

It also admits two integrals of motion since the entries of the {n + 3) x {n + 3) mass matrix M depend 
exclusively on /3o and and not on the position (a;o,2/o) of the articulated body. The coordinates xq ,yo 
are thus ignorable coordinates of the Lagrangian C and the associated linear momenta Px — dC/dio and 
Py — dC/dyo are conserved. These two symmetries reflect the fact that the system's dynamics is invariant to 
rigid translations of the whole body. Additionally, the system possesses a rotational symmetry or invariance 
to rigid rotations of the whole body. This rotational symmetry is associated with conservation of the system's 
total angular momentum h — dC/d$o + xo{d£/dyo) — yo{d£/dxo). The conservation of linear and angular 
momenta Px, Py and h can be used to eliminate the variables xq, yo and /3q from the equations of motion. 



That is to say, the {n + 3) second-order differential equations of motion obtained from (10 1 can be rewritten. 



by virtue of these three conserved quantities and a standard procedure, as (2?! + 3) first-order differential 
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equations governing the body- frame velocities UQ{t), vo{t), ^lo{t), the shape variables &{t) and their time 
derivatives ®{t). Namely, one can rewrite (10) in the form 



r) = f(r?), (11) 

where we introduced the compact notation 

r]^{uo,vo,no,&,Q}^ ■ (12) 

The function f (tj) is a nonlinear (2n + 3) vector- valued function of r] whose explicit form is omitted here for 
brevity. To get an idea of the form of these equations, note that when the articulated body is formed of a 



single ellipse with no deformation variables, that is to say, the body is rigid, ( 11 ) reduce to three first-order 
differential equations in uo,vo,no which are exactly the Kirchhoff's equations of motion for a submerged 
elliptic body; see, for example, [131 [TS]. For an articulated body with n joints, the three 'Kirchhoff-like' 
equations are coupled to 2n first-order differential equations governing the shape deformations @{t) of the 
body, n of these equations are simply stating that the time derivate of is equal to 0. 



Equation ( 11 ) admits a family of relative equilibrium solutions of the form rje = [U , 0, 0, 0,...,0]"^, 
where U is an arbitrary constant. This relative equilibrium corresponds to the motion that the articulated 
body in its elongated configuration when all its links are translating along their semi-major axis at an 
arbitrary coast speed U, with zero rotation and zero relative angles. We refer to this family of equilibria as 
the coast motion because it reminds of the coast phase in the fish burst and coast cycle. Our goal in the 
next section is to understand the role of elasticity in the passive stability of these equilibria. 



4 Stability Analysis 

We use linearization techniques to examine the stability of the equilibrium ije — [U , , ,0 , . . . ,0]^to the 
initial perturbations. We begin by presenting numerical evidence, via numerically integrating the nonlinear 



equations (111, that the system is indeed passively stable for certain parameter values. 



Before we proceed, we non-dimensionalize the equations of motion (11) by scaling the length with 
a (the length of the semi-major axis of the ellipses) and the time with a/U whereas the mass is scaled 
with pna^. All the parameters and variables can then be written in non-dimensional form accordingly. For 
example, one gets 

a=- = l, 0=-, U — — — 1 



a ' a' [/ ' * pna?U'^^ 

(13) 

_ pT:ah 7 J _ pTTab{a^ + b^)/4 _ 6(1 + b^) 

^ellipse ^ 0, e^cllipsc ;j ^ 7 • 

/97ra^ pna^ 4 

Entries in the added mass matrix can be made non-dimensional in a similar way. In this study, we assume 
I = 1.1. Note that, after this scaling, the number of independent parameters reduces to the geometric 
parameter b and the spring stiffness fc;. We want to understand the passive stability of the relative equilibrium 
fje = [ 1 1 , 0,0,... ,0]^ as a function of the dimensionless parameters b and ki. For convenience, we drop 
the tilde () notation with the understanding that all variables are non-dimensional hereafter. 

Recall that for a single rigid ellipse (see figure [ija)), this relative equilibrium is unstable for all 5 < 1. 
Due to the coupling between the translation and rotational motions in Kirchhoff's equations, an initial 
perturbation causes an elliptic rigid body to tumble when moving along its elongated direction, [Ml [15] . The 
stability of this equilibrium is more interesting when the body is deformable. 

We begin by examining the case of an articulated body made of two rigid ellipses (of non-dimensional 
minor axis b) connected via one hinge joint equipped with a torsional spring (of non-dimensional stiffness fc){^ 
The relative equilibrium in non-dimensional form is?7e = [1, 0, 0,0, 0]^. We examine the passive behavior 
of the articulated body when we impose a small initial perturbation St] to this equilibrium, for instance, 
9{0) = 0.1. Figures |3] and [4] show the response of the nonlinear equations of motion (11) for three different 



^we dropped the subscript on k because there is only one spring 
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Figure 3: Planar {xo,yo) trajectories of two-link fish for b = 0.2 and (a) k = 0.2, (b) k = 0.35 and (c) k = 0.5. Initial condition 
is Tjo = [ 1 , , , 0.1 , 0]^. Total integration time is 80. 
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Figure 4: Rotations /3o and deformations 8i versus time t of two-link fish for b = 0.2 and (a) k = 0.2, (b) k = 0.35 and (c) 
k = 0.5. Initial condition is T70 = [ 1 , , , 0.1 , 0]"^. Total integration time is 80. 
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Figure 5: Energy exchange of two-link fish for b = 0.2 and (a) k = 0.2, (b) k = 0.35 and (c) k = 0.5. Initial condition is 
T70 = [1,0, 0,0. 1,0]^. Total integration time is 80. 



sets of parameter values, h — 0.1 and k ~ 0.2, 0.35 and 0.5, subject to initial conditions rye + Sr]. Figure [5] 
shows the exchange of energy for each of these cases. Given the non-dissipative nature of the model, the total 
energy E = Tjs+Tjr + V is constant in time. But, there is exchange between the kinetic and potential energies 
of the system and consequently between the shape deformations and translational and rotational motions. 
For k — 0.35, the motion of the articulated body (figure [3[5)) remains close to the relative equilibrium and 
the body deformation Oi and rotational motion /?o (figureWfe)) remain bounded for the duration of this time 
integration. Indeed, figure[5]j6) shows that the energy exchange is just right to set the body into deformation, 
that is to say, to set 9i into oscillations, in response to the initial perturbation without causing the body 
to move far from the equilibrium solution. For the soft spring k — 0.2, the articulated body moves away 
from the equilibrium trajectory (figure [Sjja)) while undergoing large deformations 9i and large rotations /Jg 
(figure El a)) which in turn are coupled to the translational motion. This behavior indicates that the relative 
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Fi gure 6l Parameter space (6, k) of the two-link fish model. The regions for which the coast motion is linearly stable are 
highlighted in grey. 



Figure 7: Sketch of the dominant mode of oscillation in the "shark fin" stable region for the two-link model. 



equilibrium is unstable for these parameter values, figure [5ja) shows that the deformation potential energy 
V is comparable in value to the system's kinetic energy, which, by conservation of total energy, indicates that 
the coupling between the shape deformations and the rotational and translational motion causes the system 
to move away from the equilibrium. For the stiff spring k — 0.5, the articulated body also moves away from 
the equilibrium (figure [sjc)) but in a way similar to the tumbling instability observed in the case of a single 
rigid body. Indeed, the body deformation 9i is small in comparison to its rotational motion (figure |4]Jc)), 
and as in the case of a rigid body, it is the coupling between the rotational and translational motion that 
induces the instability, as evidenced from the energy exchange in figure [sjc). This is not surprising, for larger 
/c, the articulated body is more rigid and its behavior approaches that of a rigid body. Similar results are 
obtained for various values of the initial perturbations Sr] (results not shown). 

We assess the linear stability of these trajectories. We let t] = rje + Srj, where Stj is small, substitute 



r] into (11) and neglect higher order terms. The linearized equations of motion can be expressed as 

6'n = Adrj, (14) 

where the Jacobian A is a 5 x 5 matrix, and it has 5 eigenvalues Aq,Q! = 1, ... ,5. Given the conservative 
nature of the system (no energy dissipation) , the corresponding relative equilibrium is said to be marginally 
stable if all eigenvalues of the linearized system have non-positive real parts. If at least one eigenvalue has 
positive real part, the equilibrium is linearly unstable. Remember that linear instability implies nonlinear 
instability but linear stability does not necessarily indicate nonlinear stability, see, for example, |17j . For 
the two-link articulated body, we find that one of the eigenvalues is always 0, which reflects the symmetry 
associated with the fact that the equilibrium is invariant to translations in the direction of bi. The other 
eigenvalues depend on the two non-dimensional parameters of the system: the aspect ratio b of the (identical) 
ellipses and the spring stiffness k. We conduct a parametric analysis of the stability as a function of the two- 
dimensional parameter space (6, fc). Namely, we discretize a region of the parameter space from b G [0.01, 1] 



and k e [0, 2] using increments of Ab = 0.01 and Afc = 0.001, and compute the eigenvalues of (14) for each 
point of the discrete parameter space to assess the stability of the system at that point. The results are 
shown in figure [6| Two regions in the (6, k) parameter space are found for which all eigenvalues have zero real 
part, which implies that the system is marginally stable in these regions. One region corresponds to 6 > 1.68 
and the other one lies in the area 6 < 1 and resembles a "shark fin" . The existence of the first stability region 
can be understood by analogy to a single ellipse. For a single ellipse, the coast motion is stable when b > 1, 
which means a rigid body is stable when it has circular shape (b = 1) or when it is translating along its 
semi-minor axis. For the two-link case, when b ~ 1.68, the articulated body becomes more circular in shape. 
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(a) (b) 

Figure 8: Three-link fish model: (a) three-dimensional parameter space (f), fci,fc2). (b) planar cross section {b,ki = k2) 
Regions for which the coast motion is linearly stable are depicted in grey. 




(b) bending deformation 



Fi gUre 9; Sketch of the dominate modes of oscillation in the stable regions for the three-link fish model, (a) traveling- wave 
deformation mode is dominant in the small region, (b) bending deformation mode is dominant in the big region. 



granted not exactly circular, and when h > 1.68, the two- link fish model is similar to a body moving along 
its minor-axis of symmetry. The "shark fin" stability region is much more interesting. Indeed, in contrast 
to a single rigid body which is unstable for all 5 < 1, this stability region indicates that the coast motion 
can be passively stable when fa < 1 by allowing the body to deform (via the deformation variable ^i) and 
by proper choice of body elasticity k. For a given aspect ratio, if k is small (weak spring), the deformations 
are large and the system is unstable. If, on the other hand, k is too large (stiff spring), the deformations are 
small and the model is similar to that of a rigid body, which is also unstable. For an intermediate range of 
k values which provide the appropriate bending resistance to perturbations, the system is stable. 

It is informative to look at the linear modes of oscillation in the "shark fin" stable region. Remember 
that one eigenvalue is zero, say Ai = 0. The remaining four eigenvalues are two pairs of pure imaginary 
complex conjugates, say, A2.3 = ±/iiz and A4.5 — ±/Lt2i with /ii < /i2. Let ^2,3 = Ci i ^nd C4,5 = C2 ± 1772 
be the complex conjugate eigenvectors corresponding to A2,3 and A4^5, respectively, /xi and /i2 are the natural 
frequencies of oscillations and ^1, "qi and ^2, '?2 are the associated linear modes of oscillations, with the smaller 
frequency /ii dictating the dominant mode of oscillation. When h = 0.1, k = 0.35 for which the equilibrium 
is linearly stable, the smaller frequency of oscillation is /ii = 0.205 and the associated mode of oscillation is 
given by ^1 = [0 , , -0.192 , 0.257 , -0.214]^, 771 = [0 , , , -0.114 , 0]^. A sketch of this mode is depicted in 
figure [7j One can see that the two ellipses oscillate symmetrically about the hinge joint, that is to say, the 
dominant stable mode corresponds to a bending mode of deformation. 

We now ask how this stable region is affected when the articulated body has more deformation degrees 
of freedom. For concreteness, we consider an articulated body made of three identical links (of minor axis 
b) connected via two hinge joints equipped with torsional springs of stiffnesses ki and ^2- One now has two 
deformation variables 9i and 6*2 and the independent variables are r] = [uq , vq , ^0 , ^1 : ^1 , ^2 , 62]'^ ■ The 
parameter space is three-dimensional (&, fci, ^2). We follow the same procedure as above to assess the linear 
stability of the coast motion as a function of this parameter space. Interestingly, one finds two stable regions 
in the area where & < 1: a "shark fin" region similar to the one obtained for the two-link model and a 
smaller region, see figure |8] To compare this stability result with the two- link case, we plot the cross section 
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of the plane fci = k2 in figure |8|6). One can see that, in addition to the two stable regions when 6 < 1, 
a third stable region exists for b > 3.23 which is analogous to the b > 1.68 region in the two-link model. 
The "shark fin" region is also analogous to the stable region in the two-link model, though its area is larger 
than the latter. However, the small stable region does not have an analogy in the two-link model. To better 
understand this difference, we plot the dominant oscillation modes for both regions in & < 1 in figure [9j The 
dominant mode of oscillation in the "shark fin" region is similar to that observed in the two-link model. 
Namely, the two side ellipses are oscillating symmetrically relative to the middle ellipse which corresponds 
to a bending mode of deformation. Figure [9j^a) illustrates the dominant mode in the smaller region. Here, 
the two side ellipses are oscillating anti-symmetrically relative to the middle ellipse, and the deformation 
mode corresponds to that of a traveling wave. This mode cannot exist in the two-link model. An interesting 
question which remains unanswered in this work is how these stability results get affected as the number of 
hinge joints goes to infinity, that is to say, as the body becomes completely deformable. We speculate that 
similar results, with both traveling-wave and bending modes of deformations, will hold depending on the 
elastic properties of the deformable body. 



5 Conclusions 

We used a deformable body consisting of an articulated body equipped with torsional springs at its hinge 
joints as a simple model to study the passive stability of the coast motion of an underwater swimmer. The 
coast motion, that is to say, the motion of a swimmer with constant velocity along its elongated direction, 
is known to be intrinsically unstable if the body is not deformable. Through a combination of numerical 
examples and linear stability analysis, we showed that the coast motion can be passively stabilized. We 
highlight the following main points: 

(i) The surrounding fluid, even in an inviscid and irrotational context, couples the deformation of the body 
and its translational and rotational motion and allows for energy exchange between these motions in 
a way that causes initial perturbations to the coast motion to remain bounded. 

(ii) Elasticity of the body (finite k) when properly chosen can be the essential ingredient for the passive 
stabilization of an otherwise unstable coast motion. 

(iii) The effects of viscous boundary conditions and the resulting generation of vorticity from the boundary 
layer region remain an open question. One may think that viscosity and dissipation could enhance 
passive stability by damping out initial perturbations. However, it has been shown that the hydrody- 
namic forces created by vorticity shed from a bluff body may cause vortex-induced vibrations (see, for 
example, [51l28j) which may in turn have destabilizing effects. These vortex- induced vibrations may be 
especially significant where the swimmer has a bluff shape {b/a > 1.68 for the two-link and b/a > 3.23 
for the three-link). Their effects on the motion of the more fusiform-shaped swimmer {b/a < 1) is 
less clear. On the one hand, the deformation modes may add to the bluffness of the swimmer's shape. 
On the other hand, separation and shedding may be hindered by a traveling-wave deformation if the 
wave speed is faster than the swimming speed; see, for example, Borazjani & Sotiropoulos (2008), 
(2009a, b); [23l |24]. In reality, drag forces will slow down the passively swimming fish and passive 
locomotion cannot be maintained. The fish will have to actively flap in order to overcome the drag 
forces and continue swimming ~ thus the burst phase in the burst and coast cycle. During the burst 
phase, the fish could potentially correct any destabilizing forces due to vortex shedding. 

(iv) While we restricted the discussion in this paper to planar motions, the main stabilizing mechanisms 
discussed here are also present in a three-dimensional environment but for only a subset of the 3D 
perturbations, namely, for planar perturbations (see ^lOj). These passive mechanisms are not sufficient 
to stabilize out-of-plane perturbations and one needs to incorporate active control mechanisms. The 
next steps in the development of this class of models are to address three-dimensional motion planning, 
implement control strategies for both stabilization and produce complex maneuvers, such as turning 
and fleeing responses, as well as to examine multi-body synchronization (see, for example, |191I25| ) and 
its role in facilitating and/or stabilizing the coast motion. 
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